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We investigate in which cases a singular evolution with a singularity of Type IV, can be consistently 
incorporated in deformations of the R^ inflationary potential. After demonstrating the difficulties 
that the single scalar field description is confronted with, we use a general two scalar fields model 
without other matter fluids, to describe the Type IV singular evolution, with one of the two scalar 
helds being canonical. By appropriately choosing the non-canonical scalar field, we show that the 
canonical scalar held corresponds to a potential that is nearly the R^ inhation potential. If the 
Type IV singularity occurs at the end of inhation, the Universe’s dynamical evolution near inhation 
is determined ehectively by the canonical scalar held and at late-time the evolution is ehectively 
determined by the non-canonical scalar. We also discuss the evolution of the Universe in terms of 
the ehective equation of state and we show that the Type IV singularity, that occurs at the end 
of inhation, drives late-time acceleration. If however the singularity occurs at late-time, this might 
ahect the inhationary era. We also investigate which Jordan frame pure F{R) gravity corresponds 
to the nearly R^ inhation scalar potentials we found. The stability of the solutions in the two scalar 
helds case is also studied and also we investigate how Type IV singularities can be incorporated in 
certain limiting cases of R + RF gravity in the Einstein frame. Finally, we briehy discuss a physical 
appealing scenario triggered by instabilities in the dynamical system that describes the evolution of 
the scalar helds. 

PACS numbers: 04.50.Kd, 95.36.-l-x, 98.80.-k, 98.80.Cq,11.25.-w 


I. INTRODUCTION 

One of the most difficult and challenging problems in modern cosmology is the consistent explanation of singularities 
of any sort. Particularly, we need to explain if these singularities belong to the very own fabric of spacetime, or these 
are possible indicators that a quantum theory of gravity needs to take place of the classical theories that predict these 
singularities. Among singularities, there exists a classification that hierarchically determines which are catastrophic, 
or of crushing type, and which are milder. The most well known singularity of crushing type is the initial singularity, at 
which the spacetime is geodesically incomplete, meaning that the null and time-like geodesics cannot be continuously 
extended to arbitrary values of their parameters. These singularities were classified in the pioneer paper of Hawking 
and Penrose [l|, were the strong energy conditions theorems were firstly developed. On the other hand, milder 
singularities are singular spacetime points, at which the strong energy conditions are not violatech although some 
observables may blow-up at these points. These were first considered in detail in Q, but see also [j-Q for relevant 
studies, and they are known as sudden singularities. In the same line of research, finite-time cosmological sin^larities 
0 belong to the category of non-crushing types of singularities, with the only exception being the Big Rip At 

these mild singularities, the fact that some observables blow-up does not necessarily imply geodesics incompleteness. 

One of the most ’’harmless” singularity, is the Type IV singularity, which occurs if the scale factor, the effective 
energy density and the effective pressure are finite, but the higher derivatives of the Hubble rate diverge. In the 
recent studies 0113 , we extensively studied the effects of a singularity of this type, to the cosmological evolution 
of a Universe filled with one 0 or two scalar fields and particularly on certain inflationary models m- For a 
related study on this account, see 0, were similar issues are discussed in the context of sudden singularities. Using 
some standard reconstruction techniques, applicable to general scalar-tensor theories, we were able to demonstrate 
that the Type IV singularities can in some cases be consistently incorporated to the cosmological evolution of a 
Universe filled with one or two scalars and in some other cases, the existence of a Type IV singularity can prove 
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to be catastrophic, affecting explicitly the observables, as these were measured by the Planck [T2| and BICEP2 [T^ 
collaborations. The very interesting problem is then: can some (well-known) inflationary cosmologies transit to nearby 
singular inflationary cosmologies, with nearly the same values of inflationary parameters? 

The purpose of this paper is to address the issue of having a singular cosmological evolution, for a scalar potential 
corresponding to the inflation potential, or to a nearly inflation potential in the Einstein frame. In this case, 
the problem becomes more involved and as we explicitly demonstrate, the need for two scalars becomes compulsory. 
Having at hand two scalar fields, the possibilities of a successful description of the desired cosmological evolution 
becomes less involved. For relevant studies on the attributes of the two scalar field cosmological models, the reader 
is referred to [l3, [3. As we shall evince, by appropriately choosing the scalar fields, one of which is canonical 
(the inflaton) and one of these non-canonical, the Type IV singular evolution can be incorporated in the context of 
some variants of R^ inflation potential, with the inflation era being controlled by exactly this nearly R^ inflationary 
potential, while the late-time era by the non-canonical scalar field. We achieve this by abandoning completely the slow- 
roll conditions for the second scalar field, and we numerically solve the equations of motion of the non-canonical scalar 
in order to further support our claims. The model itself has some appealing attributes with regards to it’s cosmological 
implications. Particularly, as we shall show, if we assume that the singularity occurs during the inflationary era, this 
singularity affects the late-time evolution, with the latter being of nearly phantom type but infinitely close to de 
Sitter. For studies on phantom dark energy era and the possibility of crossing the phantom divide, see fl3. H^. fl8l| . 
In addition, for the possibility of having phantom inflation, see [1^ . In view of the astrophysical data that strongly 
indicate a phantom dark energy era, this could be of importance and we study this possibility in detail. In addition, 
we examine which F{R) gravity can generate the corresponding nearly R^ inflationary potentials in the Jordan frame 
and we also examine the stability of the cosmological solution, in the context of scalar-tensor theories. Moreover, 
we investigate how the Type IV singularities can be consistently incorporated to limiting cases of i? -|- R^ gravity, 
in the Einstein frame. Furthermore, we investigate which kinetic term and scalar potential govern the cosmological 
evolution in the presence of perfect matter fluids. We also provide some strong motivation for the significance of the 
Type IV singularity and finally we demonstrate that the two scalar field formalism for the observational indices can 
lead to the same results as in the single scalar field case, thus validating our approximations. 

An important remark is in order. The models we shall use are two scalar field models with one being canonical and 
the other one non-canonical. The canonical scalar part is chosen to be the R? inflation model in the Einstein frame. 
The reason for that is that the model is at 95% concordance with the current observational data. This provides clear 
motivation for using this model, and it cannot be considered as a toy model. However, the second scalar field is 
chosen in the simplest way so that a singular evolution is realized. This is a matter of choice, but for simplicity we 
chose it in the simplest way we could. In this way we shall demonstrate that, although a singular evolution cannot 
be incorporated in the single scalar field formalism, with two scalar fields this is possible. 

This paper is organized as follows: In section H, we briefly provide all the essential information for the finite time 
singularities and in section HI we analyze in detail why within the context of the single scalar field reconstruction 
method, it is a formidable task to incorporate a Type IV singularity in a i?^-like inflation potential. In section IV we 
address the problem of singular nearly R? evolution, using two scalar fields. After providing all necessary information, 
we proceed to the presentation of the model and also we find the detailed form of the nearly R? potentials in the 
Einstein frame that can incorporate Type IV singularities. Moreover we find the corresponding Jordan frame F{R) 
gravity and also we support numerically the assumptions we made during the presentation of the model. In section 
V we study the implications of our model, on the cosmological evolution, in terms of the equation of state (EoS), and 
discuss the appealing possibility that the late time era may be driven by the Type IV singularity which is assumed to 
occur at the end of the inflationary era. The stability of the cosmological solution is examined in section VI, while in 
section VH we study some limiting cases of i? -|- R^ gravity in the Einstein frame and the possibility of incorporating 
a Type IV singularity in the corresponding Einstein frame potentials. In section VHI we study the cosmological 
evolution with two scalar fields in the presence of perfect matter fluids. In section IX we demonstrate why a Type 
IV singularity is significant to study, by investigating some implications to the cosmological evolution. Particularly, 
we briefly analyze qualitatively an interesting possibility that occurs when the dynamical systems that describe the 
evolution develop instabilities, in the presence of Type IV singularities. In section X, we calculate the observational 
indices using the two scalar field formalism, and we demonstrate that the approximations we did in the previous 
sections are valid, for the choice of the parameters we made. The conclusions along with a discussion on the results 
follow in the end of the paper. 


^ Note that it is not hard to construct the model with two scalars where both scalars are of quintessence type so unifying quintessential 
inflation with Type IV singularity after it, and quintessential dark energy epoch with some of soft finite-time singularities after it. A 
review on non-singular unified evolution of quintessential inflation with dark energy era is given in 
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II. FINITE-TIME SINGULARITIES CLASSIFICATION AND CONVENTIONS 

Before we get into the problem, it is worth recalling here the classification of the finite-time singularities, according 
to Refs. [1, US]- The finite-time future cosmological singularities are classified as follows, 

• Type I (“Big Rip Singularity”) : When the cosmic time approaches t —>■ fg, the scale factor a, the effective 

energy density peff and also the effective pressure Peff diverge, that is, a —>■ oo, peff —>■ oo, and oo. For a 

detailed presentation of the Big Rip singularity, see for example Ref. Q and in addition Refs. [a@. 

• Type II (“Sudden Singularity”) d Q: When the cosmic time approaches t —>■ fg, the scale factor a and the 
effective energy density pes are finite, that is, a —>■ Ug, peff —t Ps- On the contrary, the effective pressure Pes 
diverges, that is, |peff| —t oo. 

• Type III : When the cosmic time approaches t —>■ tg, the scale factor is finite, that is, a —>■ Og, but both the 
effective pressure and the effective energy density diverge, that is, |peff| —t oo and peff —>■ oo. 

• Type IV : This type of singularity is the less harmful, with regards to geodesic incompleteness point of view. For 
a detailed account on this see d. This type of singularity occurs when, as the cosmic time approaches t —>■ tg, 
the scale factor, the effective energy density, and the effective pressure are finite, that is a —>■ Og, Peff Ps, 
|peff| —t Ps, but some higher derivatives of the Hubble parameter diverge H = a/a 

For a detailed analysis on all the finite-time cosmological singularities, the reader is referred to Ref. 0. In all the 
following considerations, the spacetime metric is assumed to be a spatially flat Friedmann-Robertson-Walker (FRW), 
of the following form, 

ds"^ =—dt"^ + ^ (dcc*)^ , (1) 

i=l,2,3 

In terms of the Hubble rate, the effective energy density Peff and the effective pressure PeS are given as follows, 

p^e = ^H\ p^g = -l-f2H + ShA . ( 2 ) 

Kj Kj \ / 

III. PROBLEM OF CONSTRUCTING SINGULAR EINSTEIN FRAME GRAVITY POTENTIAL 

WITH SINGLE SCALAR FIELD 

The purpose of this article is to investigate if it is possible to consistently incorporate a Type IV singularity to 
potentials that in the Einstein frame are nearly the inflation potential , which, as was indicated by the latest 

Planck data, can generate viable inflation [T^. Notice that the scalar-tensor potential version of the inflation were 
studied for the first time in [^ . Furthermore, the transition between non-singular nearly R^ inflation and singular 
nearly i?^-inflation remains to be the open question. The Jordan frame R? gravity in the absence of matter fluids, is 
described by the following four dimensional action, 

+ . ( 3 ) 

where with Mpi = 1.22 x 10^®GeV, and in order to have consistency with recent Planck data, we must set 

M ~ lO^^GeV. By performing the conformal transformation = /{pdjgfiy, with f{(p) = we may obtain the 

corresponding Einstein frame scalar-tensor theory, with action, 

S = I d'^xyT^ { - ^(<P) I • (4) 

with V{ip), the scalar potential which is, 


V{p) = (l - ' (5) 

Note that R? inflation from the action ([3]) is equivalent to tree-level non-minimal Higgs inflation in vacuum. 
Of course, the account of one-loop renormalization group improved corrections in non-minimal Higgs inflation [^ . 
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or even two-loop corrections completely changes the potential and breaks the aforementioned equivalence even 
in vacuum. 

In Refs. a consistent incorporation of finite-time singularities for power-law potentials was achieved by using 

one or two scalar fields. In the present work we shall attempt to do the same for the potential (IS|). In the approach we 
adopted in Ref. 00 , we made use of general scalar-tensor theories with one or two scalars, but for the potential ([5]), 
it is not easy to incorporate a Type IV singularity by using only one scalar field, as we now explicitly demonstrate. 
Therefore, the need for two scalar fields is compelling. In addition, the use of two scalar fields has very appealing 
cosmological consequences, since a unified description of singular Type IV inflation described by the potential (O, 
and of a nearly phantom dark energy era, is possible in the theoretical framework we shall use. More importantly, it 
is possible that the dark energy era is driven by the finite-time singularity it self, which is assumed to occur at the 
end of inflation. 

Before we proceed to the description of the study with two scalar fields, let us see why a single scalar-tensor theory 
fails to consistently incorporate the Type IV singularity. Consider the scalar field action given below, 

S = J I -^R - - V(</>) I . (6) 

which describes a single, non-canonical scalar field. The function a;((/>) is the kinetic function and V{(j)) is the 
corresponding scalar potential and also we assume a flat FRW background of the form ©■ Therefore, the energy 
density and the pressure are equal to, 

P = + y{(t>), P= -V{(l}). (7) 

and consequently, the scalar potential V{(j)) and the kinetic term can be written in terms of the Hubble parameter 
as follows. 


= -^H , v(</)) = 4 + h) . 

Kj \ / 

From action ®, after making the following transformation, 


( 8 ) 


P = 



(9) 


it is possible to rewrite the action ([6]) in terms of a canonical scalar field tf. Indeed, the kinetic term of the scalar 
field becomes, 




( 10 ) 


Therefore, the action ([ 5 ]) becomes. 


s = I v((p)|. (11) 

In order to study the incorporation of finite-time singularities, in Refs. 00 the scalar-reconstruction method was 
used 0, , which we now describe in brief. In the context of the scalar-reconstruction method, it is assumed that 

both the kinetic term u}{(j)) and the scalar potential F(0), are written in terms of a function f{4>), as follows, 

w(<)>) = , V{^) = 4 (3/(0)^ + /'(<^)) , (12) 

Kj Kj 

Consequently, neglecting the contribution of matter fluids, the FRW equations can be written as, 

= H = f{t). (13) 

The difficulty in incorporating finite-time singularities for the potential ([5]) is traced on the exponential form of the 
potential. In order a finite-time singularity occurs, the non-canonical scalar-tensor reconstruction function f{t) must 
certainly be of the form. 


f(t) = c{t - tsf + flit) 


( 14 ) 
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where the function fi{t) is the corresponding non-canonical scalar-tensor reconstruction function that produces the 
scalar potential m, when canonically transformed of course. In addition to (HH), there exist similar forms for the 
reconstruction function /(t), for example, 

fit) = /(t) = -h c (15) 

with c an arbitrary constant. In all the above cases, it is really difficult to incorporate these in the potential ([5]), for 
the following two reasons: 

• At first, the function f{t) must be such, so that the integral ip = ^ /'((^)d(^, can be solved explicitly in terms 

of 0 = (fip)., which is a rather formidable task for the functional forms (fTdll and (IT5|) . 

• Secondly, if someone discovers this /(</>) function, so that (j) = (f’if) is explicitly solved, then the scalar-tensor 
potential, when expressed in terms of the canonical scalar field, becomes very constrained. Indeed, suppose ipi^) 
is found explicitly, then the potential V{(j){(p)) is given by 

VW^)) = -4 (3/((/>(<^))^ + fW^))) , (16) 

Kj 

and this must be of the following form, 

Vi^i^ip)) = \m^MI (l - + V^im) (17) 

which is very difficult to achieve. 

Therefore, the need for an alternative approach is compelling, in order to incorporate Type IV or other singularities 
in the cosmological evolution. In the next section we shall make use of two scalar fields in order to achieve this. 


IV. SINGULAR NEARLY GRAVITY EVOLUTION WITH TWO SCALAR FIELDS 

The use of the two scalar field reconstruction scheme offers many more possibilities of cosmological evolution in 
comparison to the single scalar field method. In addition, possible inconsistencies that may occur in the single scalar 
field reconstruction, like for example the appearance of infinite instabilities at the cosmic time corresponding to the 
phantom- non-phantom transition, are properly amended [13 . [isL . For a recent study on this account see also M- 
We shall adopt the two scalar fields reconstruction scheme and as we shall demonstrate, apart from the consistent 
incorporation of the Type IV singularity in the theoretical framework of the nearly inflation gravity, interesting 
cosmological phenomenology is generated. Before getting into the details, we shall describe in brief the theoretical 
apparatus that we shall make extensive use of. 

The two scalar field scalar-tensor action we shall consider has the following form, 

s = J d'^xy/^ { - ^r]ix)df,xd^x - ^(<(’W) | ■ (18) 

In the above equation, the function toifi) represents the kinetic function corresponding to the scalar field (j>, and ? 7 (x) 
is the kinetic function corresponding to the other scalar field x- H is to be understood that in the case one of the 
kinetic functions oji(f>) or ri{x) is negative, then the corresponding scalar field becomes a ghost (phantom) field. For 
simplicity, we assume that the non-canonical scalar fields (j) ^^nd x depend only on the cosmic time t. For a spatially 
flat FRW metric of the form o, the FRW equation corresponding to the action (fT51) are, 

+ vix)x^ = -^H, Vif^, x) = 4 (377' + h) ■ (19) 

Kj Kj \ y 

If the generalized scalar potential ¥{(/), x) and the kinetic functions U!{(j>), rj^x) satisfy the following, 

uj{t)+r]{t) =, V(t,<) = 4 (3/(0'+/'(7)) : (20) 

Kj Kj 

then, the explicit solution of Eqs. (IT^ has the following form. 


(f = X = t, H = f{t). 


( 21 ) 
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This method materializes in brief the two scalars reconstruction scheme, a detailed account of which can be found in 
M- An attribute of this method is that in principle there is much more freedom in the choice of the kinetic functions 
u}{(f>) and ri(x)- A very convenient choice for these functions is, 

+/'(</>)"} , ^(x) = -4v«(x)^ + /'(x)^- (22) 

with the function a{x) being an arbitrary function of the scalar fields. We define a new auxiliary function /{(pjX) to 
be of the following form, 

= -y d4)u:{4>) + J dxvix)^ ■ (23) 

This function has the important property, 

(24) 

which actually fixes the arbitrary constants of integration arising in Eq. (12311 . The scalar potential can be written as 
a function of f{4>,x)j in the following way. 


mx) 


1 


^3/((/),x)^ 


^/(<(>,x) 

d(l> 


X) ^ 

dx ) ’ 


(25) 


Consequently, along with the FRW equations Eq. (HU), the following two field equations hold true, 

0 = uj{^)4> + ^uj ', 0 = 7?(x)x + ^V'{X)X^ + ^Hv{x)x + ■ (^6) 

The kinetic functions uj{(j)), rj(x) and the two-scalar potential V{4>,x), define a two-scalar field scalar-tensor theory, 
the cosmological evolution of which is given by Eq. (HIJ. In the following we shall focus to theories with a Hubble 
rate of the form, 

H{t) = h{t) + f2it)it,-tr , (27) 

with a an arbitrary parameter that actually critically determines the type of the finite-time singularity. Particularly, 
the classification of the types of finite-time singularities, for various values of the parameter a, is given in the list 
below. 


• a < — 1 corresponds to the Type I singularity. 

• — 1 < a < 0 corresponds to Type III singularity. 

• 0 < a < 1 corresponds to Type II singularity. 

• a > 1 corresponds to Type IV singularity. 

For the needs of this paper, special emphasis will be given on the Type IV singularity, which is the most harmless, 
from a geodesics incompleteness point of view. A thorough analysis for the impact of this singularity on single scalar 
field inflation was performed in [Qj. As we already noted, the function a{x) appearing in equation (|2^ . can freely be 
chosen. We shall assume that in our case it has the following form, 

Q-ix) = y [f 2 {x) {ts - x)°‘ + af 2 {x) {ts - a;)““^) - {H'{x)f . (28) 

where it is to be understood that the variable x can be either (j) or y. With this choice of a(a;), the resulting expressions 
for the scalar field kinetic functions and ri(x) are very much simplified, and these are of the following form, 

’ dix) = -4 (f^ix) {ts - x)“ + a/ 2 (x) {ts - x)“”^) , (29) 

The specific form of the function a{x) appearing in Eq. (1281) . specifies the final form of the auxiliary function /((/>, y), 
which we defined in Eq. (05)) . Particularly, this becomes. 


/(<(', X) = /i(<(') + /2(x) {ts - x)“ , 


( 30 ) 
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Thereby, the two scalar field potential has the final form, 

^ +h{x){ts-xff + f[W +fLix){ts-xT + af2{x){ts-xT~^) ■ (31) 

Having relations (l2^ and dSIl, we may easily proceed to reconstruct the deformed singular version of nearly 
inflation. Our aim is to have a Type IV singularity in the Hubble rate, while at the same time the scalar potential 
contains some deformation of the potential ([5]), at least in one of the two scalar fields. In order to achieve this, we 
assume that the Hubble rate has the following form, 

H{t) = - - -- + C 4 + C 5 (tg — t) , (32) 

C 2 + Cat 

where the parameters Ci, with i = 1, ..5, are arbitrary constant parameters. This is the simplest choice we can make 
in order to consistently accommodate the Type IV singularity in the nearly inflation cosmological model. From 
the classification given in the list above, when a > 1, the cosmological evolution develops a Type IV singularity and 
this is what we assume in the rest of this paper. In addition, for reasons that will become clear later, we assume that 
a has the following form, 


n 



(33) 


where m can be any positive integer and n is some positive even integer, chosen in such a way so that a > 1. For the 
choice (l32)) . the kinetic functions and ’qix) given in Eq. (l29t become, 


;(</)) = 


2ciC3 


K^(c2 + C3()))2 


■nix) = 


2 c 2 aits - x) 


(34) 


from which it is obvious that the scalar field x is a non-phantom scalar. The corresponding aix) function of Eq. 
now becomes. 


O^ix) = ^ clits - x) 2+2aQ,2 _ ~ " x) . 

With the choice (1551) . the function f{4>,x) becomes equal to, 

/(')', X) = -C4 + —+ csitg - x)“ , 
and consequently, the two scalar field potential V (</>, x) becomes equal to, 

vu y) = M M _£l£3_+__ 

+ ^2(c2 + c3<^)2 k2(c2 + {02 + C3</)) 

_ C 5 a{ts - x)~^+° 6c4C5(ts - x)° 6ciC5(t^ - x)° Scljts - x)^"" 

k2 K? «;2(c2 -I-C3(/)) k2 


(35) 


(36) 


(37) 


Now we perform the transformation ® to the scalar field </>, in order to transforms it, to it’s canonical scalar field 
counterpart ip. By using ([9]) and the form for a;((/)) given in Eq. ([34)l . we easily obtain the expression that relates the 
canonical scalar field p with the non-canonical scalar field (j), 


C2 + C3(j) = ev^ 

Therefore, the two scalar field action (1151) becomes. 


5 = / d^x^ I ^R - 


where the potential Vip,x) that contains the canonical scalar field </?, is equal to, 

ciC 3 e 


2 k(p k<p 

T~r/ X 3c2 3c?e cic^e 6 ciC 4 e 


+ 










CbOi{ts-x) 6c4C5(ts-x)“ 6ciC5(t5 - x)"e 3c§(t^-x)^“ 

o o ' o ' o 


(38) 


(39) 


(40) 
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and also we made use of the fact that \i a = nj ( 2 m + 1), with n =even, then a — 1 = {n — 2m — l)/( 2 m + 1), and so 
n — 2m — 1 =odd. Consequently, the following relation holds true in Eq. (123), 

(is= -(x-ts)-'+“ (41) 

The final expression of the two scalar field potential V{ip,x), when the canonical scalar field (p is taken into account, 
has similar form to some variants of the potential appearing in Eq. As we now explicitly demonstrate, the 

corresponding Jordan frame F{R) theories, are some modified versions of R^, the form of which is determined solely 
by the values of the free parameters Ci, with i = 1, ..4. 


A. Jordan frame F{R) gravity 
1. Model I: Jordan frame aR^ + R + A gravity 


Eor notational simplicity we introduce the following three constant parameters Cq, Ci and C 2 , which in terms of 
the parameters Ci, i = 1, ..4 which appear in Eq. (Uni) are defined to be, 


„ _ 3c| _ 3cf C 1 C 3 _ 6 ciC 4 

Oq — TT ^ — TT — ^5 ^2 — -o 

Using this notation, the potential (l40l) can be written in the following form. 


(42) 


Prt ^ ^ C5a{ts-x) , 6c4C5(C-x)“ , 6ciC5(C-x)“e , 3c^(C - x)^ 

yiF,x) = Vs[p>) -2-^-2--2-^-2— 


(43) 


where we have set Vs{ip) to be equal to. 


14(^5) — Co + C2e -(- Cic (44) 

In order the potential (l45ll resembles the one appearing in Eq. ®, we must set ci = |. Then, by doing so the 
potential reads, 


= Co + C2e-2v^«‘^ + Cie-v1^"‘" (45) 

As we explicitly demonstrate, the potential in Eq. (1451) can be a variant form of the potential ([S])- Before doing so, let 
us set the theoretical framework of our analysis and describe the dynamics of the scalar field y. The initial conditions 
of this scalar field can be chosen in such a way so that it starts from significantly small values of y, which during 
inflation can be negligible, in comparison to the contribution coming from the scalar field tp. In order to further 
suppress the contribution of the scalar field y during inflation, we choose C 5 <C 1 and also C 5 |ci|, with i = 1, ..4. 
Since a > 1, as the cosmic time increases, the scalar field y grows larger and we can in principle choose the parameter 
C 5 in such a way so that the contribution of the scalar field y becomes significant at cosmic times much more later 
than the ending of inflation. Moreover, we assume that the slow-roll condition does not apply for the scalar y. In 
addition to these, we choose U in Eq. (12^ . which is the cosmic time where the singularity occurs, to be exactly the 
cosmic time when inflation ends. In a later section, we shall thoroughly discuss these choices and possible alternative 
choices. In order to support our claim that the scalar field evolution goes as we just described, in the end of this 
section, we shall perform a numerical analysis to see explicitly how the dynamical evolution of the scalar field y goes. 

As a consequence of these constraints, during and at the end of the inflationary era, the potential (l43l) is approxi¬ 
mately equal to, 


Viip, x)-Co + + Cie-vT"‘^ (46) 

In addition to this, the kinetic term of the scalar y, during and at the end of the inflationary era, can be disregarded 
(for reasons we explained above), so that the action ( 1221 ) for t <ts, 


(47) 
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with Vb{}p) being defined in Eq. (1451) . So the inflationary dynamics is completely determined by the canonical scalar 
field The action (H71) . is the Einstein frame counterpart action of the Jordan frame pure F{R) gravity 






■A 


(48) 


with A = — Cl- In order to see this, we must conformally transform action (H51) . The technique is quite well 

known and for a detailed analysis on this, the reader is referred to [H, |2^ and also [s^- In addition, for a useful 
study with similar scalar potentials, see Ref. [3l|. Starting from action (l48l) . we introduce the auxiliary field A, and 
the Jordan frame action that describes a pure F(R) gravity, namely, 


can be written in the following way. 


5 = 


2k2 


S = J dx*^/^F{R) 

J {F'{A){R - A) + F{A )), 


where the F{R) function is. 




(49) 


(50) 


(51) 


It can easily be verified that by varying (ISOl) with respect to A, yields the solution A = R, a fact that validates the 
mathematical equivalence of the actions (|4^ and ([50|) . In order to find the Einstein frame scalar theory, we perform 
the following canonical transformation, 


ip = 


2fc2 


ln(F'(A)) 


(52) 


where p is the Einstein frame scalaron (or inflaton field). By making the conformal transformation of the Jordan 
frame metric, 


9^iv — 6 ^9^1/ 


(53) 


where the ”hat” denotes the Jordan frame metric, we easily obtain the following Einstein frame scalar field action. 


= J - ^ 9 ^‘'df,pd^p - V{p)^ 


A 


_ F(^) 

F'(A) F'(A)^ 


(54) 


The potential V(p) as a function of the canonical scalar field p is equal to, 
A F(A) 1 


V(p) = 


F'(A) F'(A)^ 2fc2 


g-\/2fc^/3(p^ _ I'iip p p ^g--v/2fe^/3v^ ^ 


(55) 


The function R is the solution of Eq. (j52j) . with, respect to v4, beciring in. miiid. thnt .A — H. Then, if the 

potential V{‘^) is known, the corresponding F{R) gravity may easily be obtained from Eq. ([55|) and Eq. ([5^ . Indeed 

by dividing Eq. (1551) with and taking the derivative with respect to the scalar curvature i?, we obtain the 

following relation. 


RFfi = —2 k^ 


3 d 


Vip) 


2/c2 di^ I 2(^y/2K^/3^ip 


(56) 


where Fr = Combining Eqs. (1551) and (15^ . for the potential Vg{p) of Eq. (H51) . we obtain the following 

algebraic equation. 


2CoF^ + CiFn- RFn = 0 


(57) 
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Assuming that ^ 0, we easily obtain the solution, 


F(R) 




( 58 ) 


which is nearly an R-\-aB? Jordan frame gravity. But the constant parameters Ci, i = 1, 2, 3 must be further specified 
for consistency, and as a consequence, the parameters q i = 1,..4 are further constrained. In order for the F{R) 
gravity of Eq. (l58l) to be Einstein gravity plus corrections, the coefficient of R must be set equal to one, that is, 


Cl 

2Co 


= -l, 


(59) 


As a consequence of this, in conjunction with Eq. (I42L these two relations result to ci = —C 4 . This requirement 
indicates that either Ci or C 4 must be negative. But since previously we have set Ci = 3/4, this means that C 4 = —3/4. 
Then, in order the Hubble parameter is always positive, the parameters C 2 and C 3 must be very small, that is C 2 <C 1 
and C 3 <C 1. Therefore for all times, by appropriately choosing the parameters C 2 and C 3 to be very small, the 
Hubble rate is always positive. Before we proceed to the next model, we summarize below all the constraints that the 
parameters Ci, i = 1, ..5 must satisfy, at least for this model: 


Cl = -C4 = -, C2 < 1, C3 < 1, C5 <C 1 

With regards to the later constraint, namely C 5 <C 1, specifically, C 5 must satisfy, 

3 


C 5 < 


4(c 2 + C3t) 


(60) 


(61) 


at all cosmic time values. This constraint will be of some importance in a later section. 

The situation we just studied, referred to the case that f ~ t,,, so it corresponds to the inflationary era. It is to be 
understood that as the cosmic time increases, the scalar field y, will at some point dominate during the cosmological 
evolution, while the potential of the canonical scalar field will be negligible. This will occur because the slow-roll 
approximation for the exponential potential (1451) . is violated when (f —> 00 . Therefore, the action (1391) for tg will 
take the form. 




(62) 


with Vx{<p,x) being equal to. 


H,(^,x)=^Co- 


C5a{tg-x) 6c4C5{ts-xT 3c^(ts-y)2“ 

^2 + ,^2 + ,^2 


(63) 


where we neglected the third term in the second line of Eq. (l40l) . since it becomes exponentially suppressed in 
comparison to the other y-containing terms. This late-time behavior has interesting consequences in the vacuum 
theory when matter fluids are not present. This will be the subject of a later section. 


2. Numerical analysis of the evolution of the scalar field y 


In this section we analyze numerically the evolution of the scalar y. Recall that we want to achieve an evolution for 
which the contribution of the scalar field y is negligible before, during and after the inflationary era, and dominates at 
late-time. In order to achieve this, we shall assume that the evolution of the scalar field y, does not satisfy the slow-roll 
conditions. It is worth recalling the slow-roll conditions for a canonical scalar field a. For a detailed presentation with 
regards to these issues, the reader is referred to Refs. 0|33. In the context of the slow-roll condition, the following 
constraint is assumed to hold true. 


< H(ct), 


(64) 


and in addition it assumed that this constraint is valid for an extended period of time. The constraint (1641) is known 
as the hrst slow-roll condition, and it ensures a long and finite acceleration period. In order the constraint (l64l) is 
valid for a large period of time, the following additional constraint must be imposed. 


|d-| < 


dVia) 


da 


(65) 
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which quantifies a constraint known as the second slow-roll condition. The canonical scalar field equation of motion 
in a FRW background is equal to, 


a + SHa + V'{a) = 0 (66) 

with the prime denoting differentiation with respect to a. In virtue of Eq. (j66p . the constraint (I65p can be rewritten 
in the following form, 


|d|<3il|d|. (67) 

By combining the two slow-roll conditions (j64l) and (I65L the equation of motion of a canonical scalar field (p in the 
slow-roll approximation becomes. 


tj ~ - 


1 dV{a) 


3H da ’ 

In the case that the canonical scalar field is transformed to a non-canonical scalar field x^ via the relation, 

rx 


rx _ 

CT = / Vv{x)(^x 


then, the slow-roll equation of motion for the non-canonical scalar field x, can be cast in the following way, 

5Hri{x)x + V'ix) = 0 


( 68 ) 


(69) 


(70) 


where this time, the prime indicates differentiation with respect to the non-canonical scalar field x- In addition, the 
non- slow-roll equation of motion of the scalar field x is given in Eq. (ESI), but we quote it here again for convenience, 


Vix)x + ^V'ix)x^ + ^Hr]{x)x + = 0 (71) 

2 dx 

We shall solve this equation numerically, in order to see if our argument for the evolution of the non-canonical scalar 
field Xj remains valid. We choose the parameters Ci, i = 1, ..5 in the following way, 

Cl = J, C2 = 10-4°, C3 = 10-28, C4 = -J, C5 = 10-88 (72) 

So these values of the parameters satisfy the assumptions we made previously in Eqs. dSSD and (ISTl) . In addition, we 
shall take into account that the present time in seconds is approximately tp ~ 4.25 x lO^^sec and the constant k is 
K = 8ttG = 2.0944 x 10-88GeV-8, with G being Newton’s constant. In addition, since we assumed a > 1, we chose 
a = 4/3 and also that G = lO-^^sec, the cosmic time that inflation approximately ended. There is another important 
reason behind that choice of a, since when a is of the form given in Eq. (I33L with n an even integer, and larger than 
one, the scalar field x is always a non-phantom scalar. Of course this strongly depends on the initial conditions and 
the values of the rest parameters, but this is crucially determined by the choice n =even. We shall discuss later on 
what happens in the case a > 1 and a = n/(2m 1), with n =odd integer. For the choices of the parameters we 

made in Eq. EH), it follows that the second scalar x is non-phantom for all times t. 

Using these values for the parameters, in Fig. [I] we have plotted the behavior of the Hubble rate as a function 
of the cosmic time t. As it can be seen, the Hubble rate is always positive with the choice of parameters we made. 
Having ensured that the Hubble rate takes only positive values, we proceed to solve numerically Eq. ED, in order to 
see if our assumption that the scalar field x evolves from small values during inflation, to larger values at late-time 
where it dominates the potential. We assume that the initial values at t = IQ-^^sec (approximately before inflation 
starts) for the scalar field x are x(10~^°) — 10-2° and x^(10-"’^°) — 10. As we can see in Fig. ([2]), the scalar field 
X remains small even for t ~ 10°°. In Fig. ([3]), it can be seen that after t ~ 1088sec, the scalar field x grows larger 
and larger. In addition, let us calculate the value of the scalar potential ((63l) for t = IQ-^^sec. We assumed that 
the potential of the second scalar field p dominates at that time so we explicitly verify this here. Indeed, the value 
of the scalar field x at t = IQ-^^sec is, x(10-88) = 10-2° and for this, the potential (1331) is approximately equal to 
U^((/5, x) — 10-8’^. Notice that we assumed that the scalar field (p near t = 10-8°sec takes quite large values and 
particularly those adopted in , so that all the exponentials of the scalar field p are approximately equal to one. At 
the same time, for large values of p and with the constants Ci chosen as in (ED, the scalar potential of the scalar field 
p, given in Eq. (|351). is approximately equal to Vs{p) — 18.75 x 10°, so our approximation is valid. A very important 
comment is in order. The evolution of the scalar field x behaves in the way we just described, only if it is assumed that 
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FIG. 1: The Hubble rate H{t) =+ C 4 , +C 5 {ts — t)°‘ as a function of the cosmic time with Cl = j, C2 = 10 cs = 10 
C 4 = —|, C 5 = 10“^®, tg = 10“®®sec and a = 4/3. 



FIG. 2: The scalar field’s x{i) evolution as a function of cosmic time t with ci = |, C2 = 10 C3 = 10 C 4 = — 
eg = 10“®®, tg = 10“®®sec, x(10“'‘°) 10“^°, x'(10“'‘°) 10 and a = 4/3. 


it’s evolution does not follow the slow-roll approximation. Indeed, if we assume that the field x satisfies the slow-roll 
conditions, then by solving numerically the slow-roll equation of motion for the scalar y, namely Eq. Gnu, for the 
potential (1551) . we may obtain the value of the scalar field x(t), at t = 10“^®sec, which is y ~ 10^°, so the potential is 
approximately equal to y) ^ 10®®, so our approximation is not valid. Before closing this section, let us explicitly 

see that with the choices of the parameters and the initial conditions of the scalar field y, the scalar field y is always 
non-phantom. Indeed, in Figs. |4] and [5] we have plotted the behavior of the function ? 7 (y(t)) for various time ranges. 
The first observation we need to make is that the scalar field’s y kinetic term always positive. This result crucially 
depends on the initial conditions and also to the form of a = n/(2m -I- 1), with n =even. The second observation has 
to do with the values of the kinetic term. As can be seen in all figures, the kinetic function 77(y) is negligible, even 
for t ~ 10®sec, which corresponds to a time before the nucleosynthesis and during the era that the CMB spectrum 
is fixed. In addition, by looking at Table HI we can see that rj(x) remains indeed small for a long period of time. Of 


Time 

t ~ 10 “®® 

t ~ 10 “® 

t ~ 10 ® 

t ~ 10 ®" 

»7(xW) 

1.43 X 10“®® 

3.8 X 10“®® 

9.5 X 10“® 

6.97 X 10® 


TABLE I: The scalar field’s x(t) kinetic function » 7 (x(t)) as a function of cosmic time t 

course we should bear in mind that the model we present is just a toy model, and a more concrete analysis should 
be performed, including matter fluids. Nevertheless, our assumption that the scalar y makes negligible contribution 
during the inflationary era and for a long time period after is true. In order to further support this, in Figs. [5] and 
[71 we have plotted the behavior of the scalar potential of y as a function of time. As it can be seen, the potential is 
quite small for a large period of time and also it is negligible during the inflationary era. In addition, this behavior 











13 



FIG. 3: The scalar field’s x(f) evolution as a function of cosmic time t with ci = |, C 2 = 10 cs = 10 C 4 = — 
C 5 = 10"^®, ts = 10"®®sec, x(10“'‘°) 10"^°, x'(10"‘‘°) 10 and a = 4/3. 




t (sec) t(sec) 

FIG. 4: The scalar field’s x(l) kinetic function ri{x{t)) as a function of cosmic time t with ci = |, C 2 = 10"'^°, C 3 = 10"^®, 
C 4 = -|, C 5 = 10"®®, ts = 10"®®sec, x(10"'‘°) 10"^°, x'(10"‘‘°) 10 and a = 4/3. 


can also be verified by looking at Table im where it can be seen that the potential of the scalar field x is negligible 
for a long period of time, and starts to be significant at present or much after than present time. 


Time 

t ~ 10"®'® 

t ~ 10"^ 

t ~ Itf 

t ~ 10®" 

% 

7.46 X 10"®® 

7.76 X 10"® 

1.2 X 10'® 

2.6 X 10®° 


TABLE II: The scalar field’s x(t) scalar potential Vx as a function of cosmic time t 

Before closing this section, we shall verify that indeed the value of the parameter a appearing in the Hubble rate 
of the model, given in Eq. (1321) . plays a crucial role in determining the phantom- non-phantom behavior of the scalar 
field X- Indeed, if instead of choosing a = nf{2m + 1), with n =even, we choose, n =odd, then the scalar field x is 
always a phantom scalar. This can also be seen in Fig. El where we plotted the kinetic function rj{x{t)) as a function 
of cosmic time t, for a = 5/3. 

Finally, let us note that the parameter space we used seems very much constrained, but in principle the same 
qualitative behavior appears for a wide range of the parameters. The fine-tuning imposed in the parameter space 
was necessary so that the slow-roll dynamics of the canonical scalar field ip governs the early-time dynamics and so 
the singularity does not affect the early-time dynamics, since it governs the evolution of the second scalar field x- 
Notice that the slow-roll approximation condition is very necessary for this resulting dynamical evolution. We need 
to note that when we say for a wide range of parameters we mean that C 2 for example instead of being C 2 = 10 "^*^, it 
can be chosen to be C 2 = 10"®° and C 3 , instead of being C 3 = 10"^® it can be chosen C 3 = 10"®®. So practically the 
parameters can take a continuum of values, so long as the fractions between C 2 /C 3 , C 2 /C 5 and C 3 /C 5 remain constant. 
So by fine tuning we mean exactly this and not fixing the values of the parameters to take a specific value. 
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FIG. 5: The scalar field’s xi^) kinetic function 'rj{x{t)), as a function of cosmic time t with ci = |, C 2 = 10 cs = 10 
C 4 = —|, C 5 = 10“®®, ts = 10“®®sec, x(10“"'°) ~ 10“^°, ~ 10 and a = 4/3. 




t (sec) t (sec) 


FIG. 6 : The scalar field’s x(f) scalar potential as a function of cosmic time t with ci = C 2 = 10 cs = 10 C 4 = — 
C6 = 10“®®, ts = 10“®®, x(10“'‘°) 10“^®, x'(10“‘‘°) ~ 10 and a = 4/3. 


B. Slow-roll Parameters and observational indices of the + R + K model 


Having presented the essentials of the canonical scalar field model with potential which at early times is given by 
Eq. (I45L in this section we shall extensively study the implications of this model on the slow-roll parameters and 
correspondingly on the observational indices. This analysis was also performed for a similar model in Ref. [s^, so 
the reader is referred to [s^ for details. The canonical scalar theory with potential (1^ corresponds to an i? -h B? 
gravity plus a cosmological constant in the Jordan frame, of the form given in Eq. (|4^ . By the form of the potential 
and for the values of the parameters (1721) , the potential (1451) can be further simplified to the following form, 

V/ {if) ~ Co -b (73) 

3^2 

with H = -^, since ci = —C 4 . Notice that we disregarded the term ^2 -, because for the values of the 

parameters as in GH), it is subdominant compared to the other two terms. The minimum of the potential is at (/? = 0, 
since the equation V'{ip) = 0 yields the solution = 0, for which the second derivative is equal to V" 
which is positive, so the critical point (/? = 0 is a global minimum of the potential V/((/?). Following 
assuming slow-roll evolution for the scalar field yi, the slow-roll equations for Lp during inflation read, 

with 7 = |. From Eq. (GH) it follows that quasi de Sitter acceleration is realized. The scalar field behaves as follows. 


[^, and by 






Vt 


(75) 


with to bounded at the beginning of inflation. Since during that time, the scalar field’s p values are quite large, 
the slow-roll parameters are quite small, so that the scalar field’s evolution proceeds slowly. Inflation ends when the 
slow-parameters become of order one and in order that we obtain iV = 60 e-foldings we must require that the initial 
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t (sec) 


FIG. 7: The scalar field’s x(^) scalar potential as a function of cosmic time t with ci — j, C 2 = 10 cs = 10 ca = — 
cs = 10-3®, G = 10"®®sec, x(10”‘‘°) 10“^° and x'(10”'‘°) 10 



t (sec) 


FIG. 8: The scalar field’s x(f) kinetic function ??(x(f)) a-s a function of cosmic time t with ci = |, C 2 = 10 ca = 10 
C 4 = -|, C 5 = 10“®®, G = 10"®®sec, x(10"'‘°) 10"^° and x'(10"'‘°) 10, for a = 5/3 


value of the scalar field is approximately equal to (pi ~ which is quite large (see [s^ for details). For this model, 
the primordial power spectrum is, 


A 


2 ^ 
R — 


727r2 ’ 


(76) 


while the spectral index Ug and the scalar-to-tensor ratio r as a function of the e-folding number N read. 


ris 




(77) 


Therefore, for N = 60, we have Ug — 0.9665, and also r ~ 0.0029, which are compatible with the Planck data [T^ . 
Consequently, the behavior of the model with scalar potential (|4^ is almost identical to the inflation model, with 
the only difference being the fact that the minimum of the potential is shifted in our case, since 14(0) ~ Co — 4. In 
addition, when the corresponding Jordan frame F{R) gravity is considered, the difference with the R^ inflation theory, 
is the appearance of a cosmological constant A, defined in Eq. (|15)) . As was also pointed out in [s^, the appearance 
of this cosmological constant at large curvatures needs to be explained, since it may originate from possible quantum 
effects. Regardless of that, the presence of a cosmological constant does not modify the inflationary properties of the 
canonical scalar’s field potential. 

We need to note here that the model we study in this paper, which consists of two scalar fields, and the single scalar 
field model we studied in Ref. Q, have some qualitative differences. Indeed, the present two scalar field model is 
constructed in such a way so that the singular evolution is governed by the second scalar field Xj while the inflationary 
era is governed by the canonical scalar field which has a nearly R^ potential. Therefore, in some sense the singularity 
is hidden in the x sector and does not appear or does not affect the slow-roll evolution of the first scalar field ip. This 
observation is crucial, since without the assumption of the slow-roll approximation for the canonical scalar p, the 
slow-roll parameters would not be given by the following equations. 


_ _ 1 
2k? \ v{ip) j ' K'^\V{p)j' 


(78) 


but in contrast, from the following equations, which are known to describe the so-called Hubble slow-roll parameters 
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M 


H 


SH = - 


H2' 


Vh = - 


H 

2HH ■ 


(79) 


By the analytic form of the Hubble slow-roll parameters (I79p it is conceivable that without the slow-roll assumption for 
the canonical scalar field, the dynamics of the canonical scalar field would be strongly affected by the singularity, since 
for specific values of a, the second slow-roll parameter would diverge. This could in principle have either catastrophic 
consequences for the inflationary era, if for example the singularity is chosen to occur before the end of the inflationary 
era, or lead to alternative physical phenomena, being related to the instability in the slow-roll parameters at the Type 
IV singularity. We shall not go proceed further towards this research line since this is out of the scope of the present 
paper, but a physical application of these effects is in progress. Finally, let us just comment on the possibility of 
having catastrophic consequences in the slow-roll indices, that in [^, this was the case if the singularity was chosen 
to occur before the end of the inflationary era, even in the context of the slow-roll approximation. However in our 
case, the slow-roll approximation and the fine-tuning of the parameters, protect the two-scalar field system from such 
instabilities, at least at the level of the slow-roll parameters. 


1. Models 11: Jordan frame i?"+i and R + C\R^ + C 2 {R + Ro)^^^ gravity 


The potential (1451) offers much freedom with regards to the choice of the parameters Ci, i = 1,2, 3. Following the 
lines of research of the previous section, we choose the parameters Ci to satisfy. 


Co = -Ci, C2 = - 


Co 


(80) 


Notice that ci =3/4 in all cases. The choice (1501) . would imply for the parameters Ci, i = 2, ...4 the following, 

C 4 = - 2 ci, C 3 = 6 ci (81) 

With this choice of parameters, the potential (1451) . becomes during and at the end of the inflationary era, 

+ + (82) 

where 7 an arbitrary parameter. Using Eqs. dMl), dSH) and (llll|) . we may obtain the Jordan frame F{R) gravity 
which gives rise to the potential (Hill) , which when i? ^ 7 , is approximately equal to, 


^(«) - 

which may be further approximated F{R) ~ 

By the same token, if we choose the parameters Ci, i = 1, .A to be, 


Cl 3, C 3 


Cl 


C4 = -2^1 


and by defining 7 = 9ci, then the corresponding Jordan frame gravity easily follows. 


F{R)=R+^ + 


Vs 

67 ' 36 \ 7 


-f 3 


3/2 


7 

4 


(83) 


(84) 


(85) 


C. A brief comment on Higgs inflation 


Before closing this section, we have to mention that the R^ inflation potential (1^ h as a direct correspondence to a 
certain limit of Higgs inflation model introduced by Bezrukov and Shaposhnikov [23| . So practically, the possibility 
of having a singular evolution in the context of the R^ inflation model, can be done in the Higgs inflation case too, at 
least at a certain limit. Indeed, the Jordan frame action of the Higgs inflaton model is (disregarding matter), 


S = 



^R{l + 2^n^H^H) 


( 86 ) 












with H the Higgs doublet, which is equal to, 
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(87) 


By performing a conformal transformation to the Einstein frame and by canonically normalizing the scalar field, in 
the limit where h 1, the Higgs inflaton action can be cast in the following form, 

S = j - \d,hd'^h - ^ (l - , (88) 

where the canonical scalar field h is related to the scalar field h as follows, 


h = 



(89) 


The analysis of the limiting case of the Higgs inflaton theory given in Eq. (l88l) . can in principle proceed as in the 
inflation case, which we described earlier, so we omit it for brevity. Note however, that the Higgs inflaton theory has 
the behavior of Eq. (1551) only for large field values, so it is expected that in the context of the inflation, there is 
more freedom for cosmological reconstruction model building. 


V. NON-SINGULAR DARK ENERGY ERA DRIVEN BY TYPE IV SINGULAR INFLATION 


In this section we shall study in detail the cosmological evolution of the Model I we presented in the previous section. 
Recall that the Hubble rate for this model has the form given in Eq. (l32l) . with the parameters being constrained in 
the way these appear in Eq. (IHOll . In Table I hD we have summarized the details of this model for convenience. Recall 


Hubble Rate 

Constraints on the parameters 

Singularity Type 

H{t) — 1- C4 -1- C5 {ts t) , 

Cl = C4 = 4, C2 -C 1, C3 < 1, C5 < 1, a = 

Type IV for a > 1, n =even 

Jordan Frame F{R) Gravity 

F{R) = R+£^+A 



TABLE III: Brief Description of Model I and it’s cosmological evolution 


that with the choice of parameters we made in Eq. (EH, and also with the choice of a = nj{2m + 1), with a > 1 and 
n =even integer, the scalar field x is never a phantom scalar. This however can change for a = + 1), n =odd, 

and for appropriate choice of the initial conditions, but we briefly study the implications of this case in a later section. 

It is worth recapitulate the approximations we made, because we shall make extensive use of this model in this 
section. At first, the cosmic time tg is considered to be the time at which inflation ends. At this point the two 
scalar fields action may be approximated by the action (1471) . so the inflationary era is dominated by the canonical 
scalar field g) with a scalar potential 14(‘^) given in Eq. (H51) . As the cosmic evolution proceeds and at late-time, 
the scalar field x starts to dominate the evolution, and the action becomes approximately equal to (l62l) . with the 
potential that governs the late-time evolution being I4:(7’5 x)i appearing in Eq. (1551) . This model offers a theoretical 
framework for a quite appealing evolutionary process. Particularly, as we now explicitly demonstrate, the model near 
the Type IV singularity is governed solely by the nearly R? inflation potential (1551) , with the effects of the scalar field 
X being disregarded, since these are in effect subdominant during and at the end of the inflationary era. The field x 
is assumed to have small values during that era, and we choose it’s initial conditions in such a way in order to achieve 
this. Notice that the Universe has no matter fluids present except only these two scalar fields. As the Universe evolves 
in time, then the effects of the scalar field begin to dominate, and at late-time the evolution is solely determined by 
the field x- Therefore we achieved to incorporate a Type IV singularity in a the cosmological evolution of a nearly 
inflation model, with the, apparently appealing, side effect of achieving singular inflation model during and at the 
end of inflation, and perhaps nearly phantom late-time acceleration. Particularly, the late-time acceleration maybe 
be driven directly by the Type IV singularity occurring at the end of inflation. In order to see this explicitly, we shall 
study the effective equation of state (EoS) parameter w^f /, 


WeS = - = -1 - 
P 


2H 


( 90 ) 
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which for the Hubble rate (l32)) reads, 


= —1 - 2 — 

3 (c4 + + C5(-t + ts)“) 

Near the Type IV singularity, which occurs at the end of inflation, the Type IV singularity behaves as follows, 

2(ciC3) 


( 91 ) 


WeB ~ -1 + 


3ci 


where we took into account the constraints (|60|) and also that as t —>■ tg, 

(-t + ~ 0, (-t + ts)“^0 


(92) 


(93) 


As is obvious, Eq. (1^^ describes quintessential acceleration which occurs near the Type IV singularity, and is 
controlled by the nearly inflation potential (j45|) . In conclusion we must emphasize the important result that 
inflation is not phantom but an almost de Sitter accelerating phase, with almost de Sitter meaning nearly quintessential 
acceleration. 

Correspondingly, at late-time the EoS takes the following form. 


Weff -1 - 


^ (*^4 + - C5t-b“) 


(94) 


where we used the fact that, since a is of the form (I33L with n =even, the following holds true for large cosmic times. 


{-t -f ~ (-t)““^ = 


(95) 


since a — 1 is a fraction in general, with odd numerator. Relation (1941) can describe phantom or non-phantom 
acceleration, depending on the choice of the parameters and the cosmic time value. With the choice of parameters 
we made in Eq. d72l) . the effective equation of state (iMl) . becomes, 

0.8788 10-38ti/3 4_95 10-29 

^6fT — 14” 2 ^ 

(-0.75 - 1.33 10-38H/3 + (-0.75 - 1.33 10-38^/3 + 

(_96) 

Erom this equation, it is obvious that when < ~ <p, the first term of Eq. (IMl) . is approximately equal to 1.097 x 10 ^4^ 
while the second term is 8.8 x 10“^®, so the second term dominates. It is obvious that the EoS is nearly phantom, 
since WeS — —1 — 8.8 x 10 “^®, but the phantom contribution is negligible, so practically the acceleration is almost a 
de Sitter one. The behavior of the EoS for times up to the present age of the Universe can be seen in Eig. [SI were it 
is obvious that the acceleration is practically de Sitter for all times. However, we need to stress that this behavior is 
strongly dependent on the choices of the parameters and also on the choice of a. Recall that the late-time behavior 



t (sec) t (sec) 

FIG. 9: The effective equation of state WeS as a function of cosmic time t with ci = |, C 2 = 10 “ 4 °j C 3 = 10“^®, C 4 = — 

C 5 = 10“3®, ts = 10“3®sec and a = 4/3 

for the two scalar model we just described, is governed by the scalar field x- I^i Table m we gathered our results. 
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EoS near the Type IV singularity 

nearly Quintessential Acceleration 

EoS at late-time 

Nearly Phantom Acceleration 

Scalar potential near Type IV singularity 

Va{(p), dominated by ip 

Scalar potential at late-time 

dominated by y 


TABLE IV: Summary of the cosmological evolution of Model I 


with regards to the EoS behavior. So finally, as we demonstrated, the Type IV singularity that occurs at the end of 
the inflationary era, is responsible for the late-time acceleration. Finally we have to note that this picture entirely 
changes if the singularity is assumed to occur at late-time. Then the late-time behavior is no longer controlled by 
the Type IV singularity, but the effects of the singularity might affect the inflationary era, since tg might be chosen 
to be very big. We shall present the details of this scenario in the next section in the cases that the scalar field y 
is non-phantom and phantom. Notice however that the impact of a phantom scalar on the scalar potential and also 
in the kinetic function r]{x) would completely alter the results we obtained in the previous sections, so the study of 
the phantom scalar is just included here for completeness. More emphasis was given in this section for the case of a 
non-phantom scalar and the singularity occurs at early time. 

Before going into this, we need to discuss another important issue, having to do with the value of the EoS parameter 
and the type of acceleration that the Universe experiences during inflation and at late-time. Particularly, near the 
Type IV singularity, which occurs at the end of inflation, the EoS is given in Eq. HMD, which as we stated, describes 
nearly quintessential acceleration. However, with the use of the parameter values given in Eq. dZl]), this quintessential 
acceleration is very close to de Sitter acceleration, since the second term of (IMI) . namely, , is approximately 

equal to, 2.66 x 10“^®, which is extremely small. So nearly quintessential acceleration means that the acceleration is 
very close to de Sitter expansion. The same applies for the late-time behavior, as we showed earlier, but in this case 
the acceleration is nearly phantom, but very close to de Sitter. So practically, the acceleration in this model is always 
nearly de Sitter, but note that no matter fluids apart from the two scalars, were considered in this work. 

Finally, let us discuss a vague point that may cause some confusion. In the analysis we performed we took the 
present time to be tp ~ lO^^sec. However, this estimate value should be calculated in our case by taking also into 
account the cold dark matter contribution to the total energy density. However, in our presentation we aimed in 
describing the features of the model qualitatively, at least with regards to the behavior of the EoS. The complete 
analysis should involve detailed numerical analysis, by taking into account all the observable numerical data, but such 
an analysis is beyond the scopes of this qualitative presentation. 


A. Other cases 


1. Non-phantom field case 


In the case that tg ^ tp, with tp denoting the present time, which is approximately tp ~ 4.25 x lO^^sec, the behavior 
of the resulting functional form of the EoS is altered drastically, since tg dominates over the term t. Indeed, when 
t tg, with tg denoting the time that inflation ends, the effective EoS is approximately equal to, 


Weff ~ -I - 


3 (c4 + +C5(U)“) 


(97) 


However, if we choose the values of the parameters as in Eq. and also that tg ~ 10®°sec, the EoS becomes, 


Weff ~ -I - 


1.66 (- 1 . 


- 1 . 3310 - 1 ®- 


7.5 10 ~^^ 

(10-40_|_10-2S()i! ^ 


(9.9 1041 


0.75 

10-40 + 10-28( 


which means that for t ~ 10 ®°sec, the effective equation of state becomes approximately equal to. 


Weff ~ -I -k 8.82 X I0-i°® -k 4.97 x lO’ 
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(98) 


(99) 


which describes nearly quintessential acceleration, but almost de Sitter, since the terms that drive the quintessential 
acceleration are almost negligible. In the large t regime, with tg ~ 10®°sec, when the cosmic time is nearly the present 
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time, the functional form of the EoS is given by the same relation, namely Eq. dSil), and also for the values of the 
parameters given in Eq. (EH), the EoS becomes for t ~ lO^^sec, 

Weff ~ -1 + 8.05 X f0-i°3 + 5.95 x 10-®^ (100) 

so the physical picture is the same as before, that is, nearly de Sitter but slightly quintessential acceleration. The 
same hold true if we choose ts ^ tp, with the terms that turn the EoS quintessential, being even more smaller. 

An important remark is in order. If the time that the singularity occurs is chosen to be of the order of the present 
age of the Universe, or even larger, then the contribution of the second scalar field x to the scalar potential at early 
times is not negligible anymore. A simple numerical calculation of the values of x, at t ~ 10“®®sec with the same 
values of the parameters we used before and also the same initial conditions, yields — 10“^°, and for this 

value, the potential of the scalar held x is equal to ^ 10^°, which destroys the physical picture we described earlier 
in this section. We will not go into further details towards this line of research, since the result is not so physical 
appealing. 


2. Phantom field case 

In order to study all possible cases, with regards to the behavior of the EoS, in this section we shall briehy discuss 
the case that the scalar held x is always a phantom scalar. As we already mentioned earlier, this corresponds to the 
case that a = njiflm + I), with n =odd integer and a > I. The most interesting case corresponds to ts — 10“®®sec, 
which is the time when inhation ends. Suppose for example that a = 5/3, then the EoS at late-time is approximately 
equal to, 

Wes c^-l + 1.07 X 10“®^ -k 5.95 x 10“®® (lOI) 

so the resulting picture is that, at times t ~ lO^^sec, the EoS describes nearly quintessential, but practically de Sitter 
acceleration. We have to note that this result is owing to the fact that when a = 5/3, then a — 1 = 2/3, which renders 
the term (—t + robust towards the change of the sign of (—t + ts). Finally we need to stress that the physical 

picture might completely change if we choose other values of the parameters, so these results we presented hold true 
only for the choice of the parameters given in Eq. (EH- 


VI. STABILITY ANALYSIS OF THE COSMOLOGICAL SOLUTION 


An attribute of using two scalar fields for the cosmological evolution, at least in the context of the reconstruction 
method we used, is that possible inconsistencies that may occur with the single scalar field description, do not usually 
occur in the case of two or more scalars. These instabilities were firstly observed in [T8l.[3^. and they occur at exactly 
the transition point from non-phantom to phantom and vice-versa, in the context of single scalar field scalar-tensor 
cosmology. Usually, these instabilities are infinite instabilities that is, a singularity might occur at the 

transition point. However, we must note that instabilities are not always an unwanted feature, in the presence of a 
Type IV singularity. We shall discuss this issue again in a later section. 

In this section we shall investigate the stability of the cosmological solution EH, with a Hubble rate as in (EH, 
explicitly verify our argument that the two scalar fields scalar-tensor theory guarantees the stability of the cosmological 
solutions. In order to see this, we shall rewrite the FRW equations as a dynamical system, and therefore we introduce 
the quantities and Y, which are defined in terms of the scalar fields (j) and x, in the following way, 


X,. = 




x) 

H 


( 102 ) 


Using these, the FRW equations of Eq. EH in conjunction with Eqs. EH, can be written in terms of the dynamical 
variables (II02|) in the following way, 


AX^ _ {K - l) 

AN ~ 2u}{(j))H 

^ V'ix) (^x ~ 

AN ~ 2r]ix)H 
AY _ AX^X^ (1 - y") 
AN ~ X^ + X^ 


-3{X^- Y) 

- 3 (A^ - Y) 

H X^X^ + I - Y(A,^ + Ax) 
^ m x^ + Ax 


( 103 ) 
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The reconstruction method solution of Eq. (EH, corresponds to the following values of the dynamical variables X^, 
and Y, 


x^ = i, x^ = i, r = i, (104) 

Thereby, in order to investigate the stability of the dynamical system of FRW equations (11031) . we perform linear 
perturbations of the new defined dynamical variables, of the following form. 


X^ = 1 + SX^, Xx = l + SXx, Y = 1 + SY, 
Consequently, the linear perturbations (llOSp define the following dynamical system. 


d 

div 


/ sx^\ 

( 

Ul'i'P) o 

0 

3 

= 


0 

T} (chi) 

Hvix) 

3 

\SY J 

\ 

0 

0 

-3- 


/ sx^ 
sx^ 


The matrix appearing in the dynamical system above (11061) . has the following eigenvalues, 

H 


Using the form of the Hubble rate given in Eq. 




My -- 3 - ^2 


Hvix) 

, the eigenvalues become equal to 


(105) 


(106) 


(107) 


= -3 - 

M^ = -3- 

My — —3 H- 


2c3 

(c4 + C2+C3t + C5(—t + ts)“) (C2 + Csfj)) 
— 1 + a 

(c4 + C2+C3t + {ts - X) 

ic2+cltr +C5i-t + t3)-^+^a 
(c4 + C2+C3t + + ^®)“) 


(108) 


It is clear that and My are negative for all t, at least when the parameters Ci, i = 1, ..5 are chosen as in Eq. (j60p . 
but the eigenvalue M^, can be positive when y crosses the value tg. Then it is possible that the system develops a 
saddle fixed point, which is a sort of mixed stability. This sort of instability could potentially indicate the presence 
of a new kind of physical phenomena, as we briefly discuss in a section later on. 

In conclusion, the solution (1^ . with a Hubble rate as in (15^ . is stable, for small y values, so during inflation. Note 
however that the solution m is one particular solution of the system of equations m and (1^ . In general, there 
might exist alternative solutions to this dynamical system, but the fact that the solution (l?T]) is stable makes this 
solution an attractor of the dynamical system. Therefore, any class of solutions of m and (1^ . will asymptotically 
coincide with m- 


VII. R + RP LIMITING SINGULAR DEFORMATIONS OF R + aR^ GRAVITY 

Having found a theoretical framework in order to describe the singular deformation of the gravity in the Einstein 
frame, in this section we shall discuss another interesting pure F{R) gravity, with Jordan frame action, 

5 = ^ / {R + XRP ), (109) 

where we ignored the contribution of all matter fluids and the parameter p is not necessarily an integer. By conformally 
transforming the Jordan frame action, using the method we presented in the previous section, the Einstein frame 
canonical scalar theory is. 




j -^R - - V(ip) 


( 110 ) 
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The potential V((p) appearing in the above equation is equal to, 

p 

V(p) = (111) 

with Vo being equal to, 

Vo = y (p - l)pP/(i-p)Ai/(i-P) (112) 

For an account on this model, see also [13,113 • The model described by the potential (lllll) is a deformation of the 
inflation potential, which becomes exactly the inflation potential for p = 2 . In this section we shall investigate 
certain limits of the potential (lllll) and examine how a singular evolution can be linked to the resulting scalar theory, 
with special emphasis on the Type IV singularity. The limiting case we shall be interested in is the small p limit of 
the potential V{p), which is approximately equal to. 


V (p) ~ (Kip) -i+P 


2 ^+ 2(-Vp) 3 -^+ 2 ^ (-2 + p)VoKp 2=(-i+p)3”^+2^(12 + p(-13 + 4p))VoK^p‘ 


Vo- 


-1 +p 


(-1+p)^ 


(113) 

We have to note that the small p limit is very much related to the Type IV singularity, since we are interested in 
cosmic times near the Type IV singularity. As we shall see, the approximation p —> 0, near the Type IV singularity 
is valid. We shall assume that p > 2 a choice that actually corresponds to a Type IV singular evolution, as we 
demonstrate shortly. For p > 2, and since we are interested in small p values, the most dominant term in Eq. (I113L 
is the first one, therefore the scalar potential can be approximated by. 


V(p) 


(k) -1+p 



Vop' 


(114) 


As we now explicitly demonstrate, a singular evolution with a finite-time Type IV singularity may be associated to 
the canonical scalar theory with potential given by Eq. (11141) . Indeed, consider a Hubble rate of the form, 

H(t) = fo (ts - t)“ (115) 


with a any non-integer positive number. Then, as we also discussed in a previous section, when a > 1, this corresponds 
to a Type IV singularity, which occurs at t = tg. We shall apply the reconstruction method [nils we used in section 
II, for the case of a single scalar field. Let the scalar-tensor theory be described by a non-canonical scalar field cj), 
with kinetic term u}((j}) and potential V{(j)) given in Eq. (jlll) . The by using the reconstruction method of section II, 
the potential reads, 

^(^) = ^ {3(/o (ts - - afo (ts - <(’)“”^} , (116) 

Since we are interested for t :^ts, the potential (I116|) is approximately equal to, 

{o‘fo(ts - , (117) 

In addition, the kinetic function uj((f>) corresponding to the Hubble rate (11151) is equal to, 

a;((/)) = —\afo {ts - (118) 

By transforming the non-canonical scalar to it’s canonical counterpart by using Eq. we obtain the following 
equation which relates the canonical scalar field p, with (j), 


K(a + 1) 


(ts 



(119) 


This relation validates our small p approximation, since as the cosmic time approaches the singularity, that is t ^ ts, 
then the canonical scalar tends to zero. By substituting this to the scalar potential V(cj)) given in Eq. (I116L we obtain 
the scalar potential in terms of the canonical scalar field. 


v(m) 


afo 


K (a -I- 1) 

2^2^ 


2(c-l) 

“ + 1 2(a-l) 

p a + 1 


( 120 ) 
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Thus by looking Eqs. (11201) and (lllll) . the potentials become identical if we make the following identifications, 


(k) 


P 

-1+p 




Vo 


afo / At (g + 1) ) “+i p 

[ 2^2afo J ’ P - 1 


2(a-l) 

0+1 


( 121 ) 


Thus by the second relation of Eq. (11211) . we may conclude that when p > 2, a. Type IV singularity occurs. Therefore, 
we successfully related a Type IV singular evolution to the scalar theory with canonical scalar potential given in pill) . 
What remains is to examine the stability of the solution da, which in our case is, 


</. = !, H = f{t) = foits-tr. 


( 122 ) 


The stability of this solution is the subject of the next section. 


3. Stability analysis of single scalar field ease 


As we pointed out in a previous section, when the single scalar field reconstruction method is employed, infinite 
instabilities frequently occur at the time instance that the phantom divide rceff = — 1 is crossed [l^ [l^l • We shall 
examine whether instabilities occur in the case of the solution given in Eq. p22|) . In order to study the stability of 
the solution, it is convenient to rewrite the FRW equations in terms of a dynamical system. To this end, we introduce 
the variables, and Y, which are defined as follows. 


V,. = 


Y = 


m 

H 


(123) 


Notice that, from the form of the variable V, we can realize that practically it measures the deviation from the 
reconstruction solution m- Making use of these variables and upon combining the FRW equation and the field 
equation for the scalar field cj), we obtain the following dynamical system. 


dx, /"(«(+ - l) 

dN 2f'{(j))H 


- 3 (X^ - V) , 


^ ^ n^){l-X^Y)X^ 
dN 


(124) 


where N stands for the e-folding number. We may easily verify that the solution of Eq. (US, has a direct correspon¬ 
dence to the following values of the variables X^ and Y, 


X^ = l. r = 1, 


(125) 


So practically the point (1,1) in the {X,Y) plane is the critical point of the dynamical system (I124|) . The stability of 
the dynamical system may be revealed if the system is perturbed linearly around the critical point (1,1), 


X^ = 1 + SX^, V = 1 + (5V , 

and by doing so, the dynamical system (11241) can be cast as follows. 



(126) 


(127) 


It is a textbook known fact in dynamical systems theory 
stable, if the eigenvalues of the matrix M, 


M = 


that the dynamical system (11271) can be considered 



H 


(128) 
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are negative. A simple analytic calculation of the eigenvalues yields, 

1 


M+ =- 




H{tY 


H"{t) H'jt) V 4g"(t) UH'jt) 

' H{tf H{tY 


M_ =- 


H{ty 

H"{t) H'jt) 

H{tY 


+ 3 


H"{t) H'jt) V 4g"(t) 12g'(t) 

H'{t)H{t) H{ty j H{tf H{tY 


In the case of the Hubble rate given in Eq. (11221) . the eigenvalues become in general, 
1/ - (-t + + a) + 


M.=-\ -3 


/o 


/o 


(129) 


(130) 


+ 


5(. 




+ 4(-t + ts)-2-2a(-l + q;)q! / (_i _|_ t^)-l-a(_l _|_ Q,) (_t _|_ 

/o V /o fo 


M_ = -(-3- 


{—t + ts) ^ “(—1 + a) {—t + ts) 


/o 


fo 


\ 


12(-t +t«)-i-“a 4(-t + t,)-2-2«( -l + a)Q ^ / (_^ + ^^)-i-a(_i + g) (-t + t,)-i-“a 


/o 


/o^ 


fo 


fo 


Although by appropriately choosing a, and also by assuming that t is away from the singular point, the above 
eigenvalues can be negative. But as the cosmic time approaches tg, the eigenvalues have an infinite instability, with 
the eigenvalues being positive in some cases. This makes the system unstable at the transition point. Note however 
that instabilities may serve as indicators that the evolutionary process may change at the time these occur, as we 
briefly discuss in the next section. 

Before we close this section, for completeness we shall briefly provide a description with two scalar fields, one of 
which is a non-canonical field. We presented the essentials of the two scalar field method in a previous section, so the 
action is given by Eq. ®, the kinetic functions by (|2^ . while the scalar potential is given by Eq. (EH). By choosing 
the auxiliary function a{x) as follows, 


a{x) = fo{-x + ts) 

we get the kinetic functions of the scalar field (j) to be, 

2/oa(tg - 0)-!+“ (1 + V2) 

W = -2-^^ 

while the kinetic function of the scalar field x reads, 


vix) 


2V2/oa(tg - x)~^+" 


The corresponding scalar potential V{(j)^x) is given by, 

. _ yboafts - + boa{ts - 

3 (6o(l + a){ts - (())“ + 2y%a({ts - 4>)^ + (ts - x )^)) 
^ (1+0)2^2 


(131) 


(132) 


(133) 


( 134 ) 
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VIII. INCLUSION OF MATTER FLUIDS IN THE COSMOLOGICAL EVOLUTION 

Having studied the cosmological evolution of Eq. (1321) in the absence of matter fluids, in this section we shall 
include the contribution of perfect matter fluids. Particularly, we shall include the effects of pressure-less matter, but 
in principle other matter fluids can be considered too. We start off with the description of the equations of motion in 
the presence of matter fluids, with constant equation of state Wm, which take the following form, 

= Y ^ + \v{x)x^ + V{(j), x)^ , (135) 

2 

^ = -y (pm +Pm +w('(>)<^^ , 

where pm and Pm stand for the effective energy density and effective pressure of the matter fluids with equation of 
state Pm = WmPm- The effective energy density and pressure satisfy the conservation law. 

Pm + iH {pm +Pm) , (136) 

from which we get pm = Now consider the following solution to the FRW equations (I135|) . 

c^ = X = t, H = f{t). (137) 

Then then kinetic terms uj{(j)) and 77(x), satisfy, 

w(0 + Xit) = -\nt) + {wm + 1) , (138) 

where F{t) = The last term follows from the fact that pm = and also by combining the first 

relation of Eq. (11351) and H = a/a. So in general, the kinetic terms that generate the cosmological evolution (I137L 
can be chosen as. 


+ (139) 

^(X) = - y (v'^^W + TW) + + . 

where ai{x) is an arbitrary function. Indeed, it can be easily verified that the above two kinetic terms satisfy Eq. 
(I138L when cj) = x = t. Also we define the function, 

/(<(>, X) =-y (^y w(<?i)d^!)-f J v{x)^X- J j /i(x)dx^ , (140) 

with /i(x) = (a:)_ function f{(l),x) satisfies = f{t), as it can be easily verified since, 

J (w(0 + V{t) - hit) - hit)) ^t = J fit)dt = f{t ), (141) 

which holds true, due to Eq. (I138|) . So the potential that generates the cosmological evolution (11371) . can be chosen 
as follows, 


_ {Wm - 1) ^ 

We chose the potential as in Eq. (11421) . so that at 0 = x = ^ it satisfies. 


v{t,t) = ^ {f{tr+fit)) 

Ki 


jWm l) ^^^_3(i+. 




(142) 


2 


(143) 
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It is easy to prove that the potential V (</>, x) of Eq. (11421) . since at i/> = x = t, the following relations hold true: 

3/(t,t) = 3/(t)2, (144) 

df 


d(j) 

dx 


= uj{t) - 




= x{t)- flit) 


and consequently at y = ((i = t, we have, 


dij) 


4>=t dx 


x=t 


= r]{t) + uiit) - 2/i(f) =-, 


(145) 


where we also used Eq. (jl38ll . Let us find the corresponding kinetic terms w(<()), r]{x), the function / and the potentia 
V{(j),x), for the cosmological evolution of Eq. (15^ . By combining Eqs. (15^ and (11391) . we get the kinetic terms 
and »7(x)j which are, 


, , 1 -3(l+-uJm) ^ +C4b+|k In(e2+C30)') 

Ujicj)) = -e V 3 y 

1 + -|_c 4;^+£1 ln(c2+C3X)') 

Vix) = 2 ® ^ ^ ) 

while the potential of Eq. (11421) . becomes equal to. 


Eb(l + Wra) + 

Fo(l + Wra) - 


2 ciC 3 


K^(c2 + C3(())2 
2c2Q:(ts - 


(146) 


— _lg-3(l+’"m)(c4b+C5(t„-<^')“(-i^ + Y^) + || In(c2+C30))^^ 

_ lg-3(l+u'm)(c4X+C5(i«-x)“(-T^ + T^) + |j In(c2+C3x))^^ 

lg-3(l + U)m)(c40+C5(t4-0)“(- 4^ + j^) + |i In(c2+C30))^^^ 


(147) 


lg-3(l+u'm)(c4X+C5(ts-x)“(-T^ + Tfe) + ly ln(c2+C3X)) 

4 

^ 3cj _ C1C3 


FoWm 


6 ciC 4 




ft: 2 (c 2 + 030)2 k 2 (c 2 +C 30)2 ^ 2 (^ 2 + 030 ) 


CbOlits-x) , 6c4C5(ts-x)“ , 6ciC5(f3-x)“ , 3c^(ts-x) 


2a 


k 2 «:^(c 2 +C 30 ) k 2 


We need to note we chose the function ai(x) as in Eq. (ESI- Eor an application of the method we used to include 
matter fluids, but in the case of a single scalar field, the reader is referred to Ref. [^ . 


IX. SIGNIFICANCE OF THE TYPE IV SINGULARITY 


In the previous sections we investigated how a Type IV singularity can be connected to the cosmological evolution 
of the E? inflation model. As we discussed this can be done only with the use of a second scalar field, which is 
dynamically insignificant at early-time. However, we did not discuss our motivation for studying such a singular 
evolution. In principle it could be claimed that a Type IV singularity is not a real singularity, but now we provide 
some arguments to support our study and the significance of the Type IV singularity. Firstly, it is not a crushing type 
singularity so the Universe can smoothly pass such a singular point and continue undisturbed its evolution. Secondly 
and more importantly, the Type IV singularity can affect the observational indices of inflation in a dramatic way, 
if the slow-roll condition is abandoned. As we will demonstrate, the effect of a Type IV singularity is to cause an 
instability to the dynamical system, which can indicate that the attractor solution that drives the early-time evolution 
is unstable, and in turn this could be an indicator that inflation ends and graceful exit is achieved. In the rest of this 
section we thoroughly discuss this issue, but for a similar situation consult Ref. [s^. 

We shall demonstrate that the Hubble second slow-roll parameter rjH develops an instability at the point where a 
Type IV singularity occurs. Recall that the second Hubble slow-roll parameter is equal to [3^ : 


Vh = 


H 

2HH 


(148) 
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and in inflationary dynamics, the second Hubble slow-roll parameter measures how much does inflation lasts, while 
the first Hubble slow-roll index measures if inflation occurs in the first place. Note that we abandon the slow-roll 
condition for all the scalar fields of our model. Let us calculate the second Hubble slow-roll index for the Hubble rate 
appearing in Eq. (15^ . assuming that a = 4/3 or more generally that 1 < a < 2, as in the previous sections. For the 
Hubble rate of Eq. (1321) . the second Hubble rate (11481) becomes. 




m = 


(c4 + + C5(-t + tsY 


(C2-|-C3i) 


T — C5(—t-f ts) 


(149) 


By looking Eq. (I149L we can easily observe that the term ^ (—f -|- appearing in the numerator of the 

fraction, becomes divergent at t = ts, for the values of a we assumed. This clearly indicates an infinite instability of 
the inflationary dynamics. In order to make things more clear, let us recall the significance of the inflationary indices 
and their interpretation. We follow the analysis of Ref. 37|. As we demonstrate, this instability could be viewed as 


an obstruction in the inflationary evolution and hence can act as a mechanism for graceful exit from inflation, or at 
lest as an indicator of exit from inflation. However, this mechanism is different in spirit from other mechanism for 
graceful exit, nevertheless the instability occurs, so we should analyze it’s significance. 

A clear interpretation of the slow-roll indices was provide in Ref. [3^ , where the slow-roll expansion was used. The 
inflationary dynamics is determined by the slow-roll expansion, and particularly the slow-roll expansion indicates if 
an inflationary attractor is the correct inflationary solution. It is a perturbation expansion, the first terms of which 
are the usual slow-roll indices. Note that the slow-roll expansion is a more strong and restrictive physical description 
of inflation, in comparison to the slow-roll approximation [3^ . Let us briefly describe the slow-roll expansion by using 
a general example. Consider an inflationary solution described by a Hubble evolution with ip being a canonical 

scalar field. The Hubble slow-roll expansion determines an inflationary solution which is an asymptotic attractor of 
all the inflationary solutions obtained by the usual potential slow-roll approximation. The FRW can be written as, 


n , , OTTK , . / i 

H i^) = —- gC// 
where en is the first Hubble slow-roll parameter. 


(150) 


= (151) 

but in Eq. (11501) is expressed as a function of the scalar field (p. The Hubble slow-roll expansion is obtained if the 
FRW equation is expanded in a perturbation series, 

+ €v- ^4- + (152) 

+ f 4.. + jevxt - f 4-).^ + + o.) , 

where terms up to fourth order are kept, and the parameters ev,77y,^v are given in terms of the Hubble slow-roll 
parameters tjh and ch as follows, 

W = (3 - Eff(Seff + -Vh- Cff), (153) 

= (3 — e_f/) ^ (^27eH'r]H + ~ ~ ~ ■ 

In addition, and an are defined in terms of en and tjh below. 


= ^hVh - y - Vh) - y■ (154) 

In the above equation, the prime indicates differentiation with respect to the scalar field ip. It is conceivable that 
the Hubble slow-roll expansion we just described is a better approximation to the final inflationary attractor H{t), in 
comparison to the standard slow-roll approximation, which uses only the parameters e and rjv at lowest order. Note 
that the potential slow-roll parameters correspond to the ones defined in Eq. (1751) . in terms of the potential of the 
canonical scalar field. In addition, a useful relation that we will make use in the sections to follow, is that at lowest 
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order, the Hubble slow-roll parameters eHiTlH and the potential slow-roll parameters ev,? 7 y for a canonical scalar 
field are related as follows, 


eff = ey, r]H = r]v - e. 


(155) 


Coming back to the Hubble slow-roll expansion interpretation, as was also discussed by the authors of 37[ , when the 
slow-roll perturbative expansion breaks down, the inflationary solution ceases to be the final attractor of the theory, 
thus inflation ends. This breakdown of the perturbative expansion can occur when the slow-roll parameters take large 
values, or if a singularity occurs. The latter is the case for a Type IV singularity. Hence our motivation for studying 
the Type IV singularity is exactly this non-trivial feature of the inflationary dynamics. 


X. OBSERVATIONAL INDICES WITH THE DOUBLE SCALAR FORMALISM 


In this paper we made a crucial assumption, with regards to the early-time dynamics of the scalar fields Lp and x, 
which was that the scalar field x does not make a significant contribution at early-time. Thus the contribution of the 
scalar field, in the potential and also its kinetic term, can be safely disregarded. This result was supported numerically 
in the previous sections, however, we need to further support this by using the multi-scalar field formalism for the 
calculation of the slow-roll and observational indices. We need to stress that when both scalar fields are taken into 
account, the space of trajectories is much more large in comparison to the single scalar field case. With the present 
investigation we aim to demonstrate that when the parameters are of the same order as the ones given in Eq. ([721). 
then the observational indices and slow-roll parameters lead to the same result as in the single scalar field case. We 
shall use two approaches with regards to the slow-roll parameters, the one developed in Ref. 38|, which is related 


to the slow-roll condition, and the second one is related to the so-called Hubble slow-roll parameters, developed in 
Ref. [S^. Note that the first slow-roll parameter e coincides in both the methods we shall use. However the second 
slow-roll parameter is different in the two approaches. For more detail consult Ref. 371. 

We use the formalism and notation of Ref. [Sa . We consider the multi-scalar field action, 


S = 


d" 


xy-g 


R 




(156) 


with I,J = 1,2. In our case the scalar fields 4>^ are (j)^ = ip and (j)^ = x- Moreover, the metric of the scalar field 


configuration space Gij{(j)^), depends on the scalar fields. In our case, 
matrix representation of the metric for the two scalar field action of Eq. 


the field space is two dimensional, and the 
equals to, 


G = 


1 0 

Q 2c2a(x-ta)~’'^° 


(157) 


The FRW equations for the action (11561) are. 




■ ^(y) 


(158) 


I 


H = - 

□y + - g^^Vk = 0, 


with Vk = dV/d(j)^, and also with we denote the Christoffel symbols corresponding to the two 

dimensional scalar field configuration space, with the metric being the one given in Eq. (11571) . For this metric, the 
only non-zero Christoffel symbol is, 


rx 

XX 


— 1 -(- o 

2(—ts + x) 


(159) 


For the sake of notational simplicity, we introduce the scalar field a and the vector field which are defined to 
be. 


d = V Gij<j3^4i-’, 

a 


(160) 
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In terms of & and , the field equations (11581) are written as follows, 


a + 3H& + = 0 

where Va = ^^Vj. In our case, the potential V of Eqs. (11581) and (11611) is identified with the potential of Eq. (1401) . 
Since a will be extensively used in the following, it is worth computing it, and for the metric of Eq. (jl57L it reads. 


a = 



2c2a(-t« + x)-i+“ 


-X 


(162) 


Up to this point, the potential slow-roll approach with two scalars and the Hubble slow-roll approach is the same. 
However, we differentiate our analysis here, since the second slow-roll parameter is different when calculated in the 
context of the two aforementioned approaches. We start of with the potential slow-roll parameters. 


4- Potential Slow-Roll Parameters Approach 


We calculated the slow-roll indices in the potential slow-roll approach. We use this name, since as was shown in 
[S^, the indices are related to the scalar potential directly. However, we shall use the slow-roll limit for the canonical 
scalar field (p, and also we shall take into account the values of the parameters as in Eq. dZl]). As we now demonstrate, 
the observational indices are approximately equal to the ones corresponding to the single scalar field case. 

In both cases, the first slow-roll parameter e, is equal to. 


H _ 


3o-^ 


a^ + V 

We can easily calculate this for the potential of Eq. (go]), so the parameter e reads. 


e = 3- 




-X 




(163) 


(164) 




For the values of the parameters appearing in Eq. dZl]), we can simplify the parameter e, since the following hold 
true. 




2c2a{-ts + x) 


— 1+a 


• 2 *2 

-X -p , 


(165) 


V{p,x) v{p) 

with V{p) given in Eq. (HSl) . Therefore, it simply follows that e equals to, 


e = 3 




= 3^ 


V 

2 




(166) 


Recalling that we assumed that the canonical scalar field p satisfies the slow-roll conditions, then as was shown in 
( 3 ^ . the parameter e equals to. 


e = ch = 5 


^ j_ ( rip) \ ^ 
l^ + V{p) 2 k 2 \ V{p) J ■ 


!£- 

2 


(167) 


Note that this results holds true at first order, as we briefly explained in the previous section. Now we proceed to the 
second slow-roll parameter rj. We denoted this as 77 in order to discriminate this from the Hubble slow-roll parameter 
rjH, which we compute in the next subsection. Since the parameter e takes small values, as was shown in [s^, the 
potential slow-roll index 77 is approximately equal to, 


77 ~ e 


fj 

lid 




(168) 
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We can easily calculate a from Eq. (I160|) and it reads, 




so in our case it is equal to, 

4 c 2(-1 + a)a{-ts + 

2 _ 2c2a(-ta+x)~^+°X^ 


+ 


2JGijcpicpJ 


2ipip 


k^ 2G tpp- — 


2 _ 2c2a(-ts+x)~^+°‘X^ 


2G^^- 


Ac2a{-ts + x) ^+“XX 


2 _ 2e2a(-ta+x)~^+°X^ 


IfP' — 


(169) 


(170) 


As in the previous steps, for the choice of the parameters given in Eq. dll, the parameter ij can be approximated by 
a ~ (p. So finally the second slow-roll index r\v reads. 


?7~ e- 


Hip 


(171) 


But as was shown in [s^, the second term of Eq. (I171|) . is the second Hubble slow-roll index for the canonical scalar 
field, that is, rjH = and therefore the we get, 


r] = e + r]H. (172) 

As was explained in the previous section and also was shown in [s^ , the expression e + r]H is equal to the potential 
slow-roll parameter rjy corresponding to the single scalar field ip (see Eq. (11551) 1. Therefore we proved that with the 
choice of the parameters made in Eq. dZH), the two scalar field potential slow-roll index rj is actually equal to the 
single scalar field slow-roll parameter rjv , which we defined in the previous section. This result validates our original 
claim that the two scalar theory is dynamically equivalent to a single scalar theory. In practise, we chose a trajectory 
for which the potential and kinetic term of the field x could be disregarded, which is equivalent to a single scalar field 
trajectory. Therefore, at the observational indices level, the single scalar field trajectory and the two scalar theory 
with the parameters chosen as in Eq. dZH). However potential differences between the two approaches can be revealed 
because the small contribution of the scalar field x can generate non-Gaussianities. This study is quite interesting 
but exceeds the purposes of this article. Before closing this subsection we need to note that the spectral index of 
primordial curvature perturbations and the scalar-to-tensor ratio, corresponding to the slow-roll indices we calculated 
in this section, are equal to. 


= 1 — 6e -I- 2?7, r = 16e , 

For the Hubble slow-roll parameters this relation is different, as was shown in [s^- 


(173) 


5. Hubble Slow-Roll Parameters Approach 

In this section, we adopt a different approach in comparison to the one we adopted in the previous section. Specif¬ 
ically, the main difference will be in the definition of the second slow-roll parameter rj, since the definition of the 
first slow-roll parameter coincides with the one we used previously. Particularly, we shall use the Hubble slow-roll 
parameter tjh, which is equal to [3^ : 


Vh = - 


H 

2HH ' 


(174) 


We can easily write the Hubble slow-roll parameter r]H in terms of the a held we introduced in the previous section, 
and by using the FRW (I161L we get. 


H = -K^aa. (175) 

Then by combining Eqs. (I175|) . ()161l) and ()169|) . the second Hubble slow-roll parameter r]H can be written as follows, 


Vh = 


&H 


(176) 
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By using Eqs. (I162|) and (11701) . the parameter t]h can be written as follows, 


Vh = 


-Ac 2 a{-ts + x) + 2(p<f 


2H 




2c2a(-ts+x) 

-- 


- 4c2(-1 + a)a{-ts + x)"2+“x 


(177) 


and by using the values of the parameters given in Eq. dZl]), the parameter tjh can be simplified as follows, 


rjH ^ - 


Hip 


(178) 


By looking Eq. (I178p , we realize that this is nothing else but the Hubble slow-roll parameter corresponding to a single 
canonical scalar field ip (see previous section). Therefore in this case too, the two scalar-field formalism ends up to 
the single scalar field description, if the values of the parameters are chosen as in Eq. ([72]). Note that in the present 
case, the spectral index of primordial curvature perturbations Us is equal to [s^, 


rig = 1 — 4e -I- 2r]H ■ (179) 

Also we need to stress that if we substitute the first order in the slow-roll expansion result uh = + ??, then Eq. 

(I179|l becomes identical to Eq. (117311 . 

As we already stated, this result is strictly dependent on the fine-tuned choice of the parameters (17^ . Therefore, if 
we can safely neglect the kinetic term and the contribution of the scalar field x to the potential, then our results are 
valid. However, some differences between the two scalar field theory and the single scalar field theory can be observed 
when non-Gaussianities are considered. We defer this task to a future publication. 


XI. DISCUSSION AND CONCLUSIONS 

In the previous section we mentioned that within the theoretical framework of general scalar-tensor theories, if the 
FRW equations are combined with the scalar field equation of motion, these can be written in terms of a dynamical 
system. The solution of the reconstruction method, is a solution to this dynamical system, so the stability of this 
solution may reveal if the reconstruction solution is stable towards linear perturbations. Here, we shall discuss the 
possibility of having an instability in the dynamical system and we shall offer another perspective of having an 
instability in a theory containing a finite time Type IV singularity. In principle, the existence of instabilities may be 
considered as an unwanted feature, but the existence of the instabilities may possibly indicate the point at which a 
new physical phenomenon takes place. This is true, since at the instability point, the classical trajectory that the 
dynamical system followed is unstable, and therefore the system becomes strongly unstable. The important thing 
about the Type IV singularities is that the dynamical system of the scalar fields becomes unstable at exactly the time 
at which the singularity occurs. Therefore, if we appropriately choose this time, it is possible to explain indirectly 
why physical phenomena described by the dynamical system, at the point of the singularity need a new physical 
description yet to be found. An example on this account could be for example the stopping of inflation in certain 
scalar theories of inflation. We hope to address this issue soon and work is in progress. 

In conclusion, in this paper we constructed a theoretical framework consisting of two scalar fields, in order to 
support a singular evolution of the Universe with a nearly inflation potential. The singularity we incorporated 
in the evolution is a Type IV singularity, which we assumed it occurs at the end of the inflationary era. The model 
we used consisted of two scalar fields, one canonical and one non-canonical, with the canonical scalar being the one 
describing the nearly inflation potential. Near the singularity, and by appropriately choosing the parameters, the 
canonical scalar field dominates the evolution of the Universe, since the non-canonical scalar field is chosen in such 
a way so that it’s contribution at early and intermediate times is negligible, in reference to the contribution of the 
canonical scalar. In addition, we assumed that the non-canonical scalar does not satisfy the slow-roll condition, and 
therefore it’s contribution becomes significant and it dominates the Universe’s evolution at late-time. As we showed, 
in some cases, the scalar model we used has some qualitatively appealing attributes, with regards to the cosmological 
evolution. Particularly, it is possible that the Type IV singularity that occurs at the end of inflation has a direct 
impact on the late-time dark energy era. Specifically, the late-time acceleration is driven by the Type IV singularity 
and moreover the dark energy era is a nearly phantom era, infinitely close however to a de Sitter expansion. 

Having found which nearly R^ inflation potentials can incorporate the Type IV singularity in their theoretical 
framework, we found the corresponding Jordan frame pure F{R) gravity and also worked out the stability of the 
cosmological solution of the reconstruction method we used. In addition, we studied certain limiting cases ot R + RP 
theories in the Einstein frame and investigated how these theories can incorporate a Type IV singularity. 
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In fact, we demonstrated that singular inflation is quite possible, while transitions between singular and its com¬ 
panion non-singular inflation is not easy to realize. In the absence of a complete quantum theory of gravity, which 
simultaneously incorporate the singularities in a consistent, complete and unified theoretical framework, we believe 
our work serves as another step towards the understanding of the nature of these singularities and more importantly, 
of their impact on the Universe’s current evolution. In view of the fact that what we actually see now is the remnants 
of the primordial quantum theory of gravity, without the quantum phenomena playing an important role, the classical 
cosmological finite time singularities are useful tools to explore the quantum phenomena within a classical theoretical 
framework. More importantly, since Type IV singularities are not-catastrophic singularities, like the initial singularity 
or the Big Rip finite time singularity, these may serve as a doorway to the complete noesis of the quantum nature 
of cosmological phenomena. In fact, in Loop Quantum Cosmology [s^, the initial singularity does not occur, but 
finite time singularities may occur. The Type IV singularities do not cause geodesics incompleteness of spacetime, 
so in principle they are harmless, but however the complete nature of these has to be understood, since these may 
have direct impact on observational indices It is possible that these singularities may be responsible for ending 
acceleration or deceleration eras, through instabilities they cause. The indicators of this effect may be the instabilities 
that often occur in the theories described by one scalar or even two scalar fields. So practically an instability that 
occurs at the time when the sin gula rity occurs may practically signal the change of the evolutionary mechanism at 
that point, see for example Ref. (^ . 

Finally, a last remark for the model we used in this paper. Basically, it is the B? inflation model in the Einstein 
frame plus an arbitrary but the simplest deformation, that makes the evolution singular. The models we studied 
cannot be considered as toy models, since these are realistic models, with regards to the canonical scalar field </?, and 
the X part can be chosen in a simple way so that singular evolution is realized. So practically the (p was chosen on 
the basis that the inflation model is in agreement with current observational data [T^ . As we demonstrated, at 
the level of observational indices, the scalar field x does not contribute significantly at early times, while at late-times 
it dominates the evolution. A quite interesting task to investigate however, is to investigate the effects of the second 
scalar fields on the non-Gaussianities of the model. Particularly it is possible that the second scalar could generate 
non-Gaussianities at early-time, although it makes minor contribution to the early-time evolution. We hope to address 
this issue soon. 
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